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Theory for Some Asynchronous 
Time-Division Switches 

By J. E. MAZO 
(Manuscript rweived September 24. 1970) 

A two-wire active asr/nchronous lime-division switch kas recently 
been proposed by J. 0. Dimmick, T. G. Lewis, and J. F. O'Neill. The 
interrupted energy transfer from one filter to another is accomplished 
asynchronously in order that more efficient use may be made of pro- 
cessing time of talker pairs on the switching bus. After first formulat- 
ing, in q concise mathematical fashion, the effect of posting a signal 
through such a randomly time-varying circuit, we foais attention on 
optimizing an important filter response function. Typically, two per- 
cent rms jitter in the transfer times yields an output S/N of SO dli 
independent of signal spectrum. The timing stabilization required to 
obtain small jitter is also discussed and an exact solution for ex- 
ponential processing time is obtained. This latter result may be put 
to good use in studying the efficiency of the switch. Conservatively, 
asynchronous operation should increase traffic capacity threefold. 

Finally, a speech wave which passes through a sampte-and-hold 
circuit with random sampling times is considered upon being recon- 
structed with a fixed filter. This is a model for a four-wire active 
asynchronous sicitch. and results are compared with the two-wire 
situation. 

I, INTRODUCTION 

Voice switching systems have most commonly been based on con- 
trolling electromechanical switches which select and hold a spatially 
distinct path for each conversation. The technology used to implement 
such a space-division network (crossbar or ferreed switches! usually 
results, in practice, of an individual path having much larger band- 
width than is required for faithful transmission of the signal. The 
space and cost of these switches makes other solutions desirable for 
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many applications. One line of attack has been to keep the space-divi- 
sion concept, but replace the electromechanical relays with semicon- 
ductor switches. These techniques, however, still suffer from the 
hard-to-grow nature of multistage space division networks. A more 
promising solution seems to be tb place all the conversations on a 
wideband bus using time-division techniques. In fact, the 101 Elec- 
tronic Switching System (ESS) is such a time-division switch. This 
system uses resonant energy transfer to "move" periodically measured 
samples of a speech wave from an incoming line to an outgoing one. In 
the same vein, an asynchronous time-division switch has recently 
been proposed by J. O. Dimmick, T. O. Lewis, and J, P. O'Neill. 1 This 
switching arrangement makes use of active energy transfer between 
filters rather than resonant energy transfer and allows a variable time 
slot for transferring each speech sample through the switch. The asyn- 
chronous nature of this switch allows a more efficient use of processing 
time than is possible to achieve synchronously. However, a consequence 
of this virtue is a periodic sampling of the input waves. While the 
synchronous switch can make use of the usual sampling theory to 
guarantee faithful reproduction, the asynchronous switch cannot. 
Further, the modifications of the theory which are required to discuss 
the proposed asynchronous switch are not simple, for the random 
sampling causes some feedback energy which is later retransmitted, 
further adding to the output noise. Our immediate purpose will be, 
then, to present theoretical work relevant to this problem. We discugs 
the quality of transmission [measured by the output signal-to-noise 
ratio (S/N) | as a function of jitter in the sample values and as a 
functional of a certain filter response function upon which the feed- 
back energy depends. The optimum function is found. Also a tech- 
nique suggested in Rei 1 for keeping the jitter small is discussed 
theoretically, and the combined question of how jitter, quality of 
transmission, and increased efficiency are related is answered. 

Section X summarizes our conclusions and, barring some terminology 
introduced in the text, may be read next. 

II. MATHEMATICAL MODEL 

Consider the diagram in Fig. 1 which represents a talker and listener 
on a switching bus. There will be many such pairs on a particular bus, 
but we need now concentrate on only one. The way the switch works 
is that, at approximately periodic instants of time, the two identical 
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Fin, 1— Model for a "talker"— "listener" pair on an asynchronous switching bus. 

capacitors which are shown in Fig. 1 are interchanged* and energy 
transfer is effects The following assumptions are made: 

(i) The input signal x(t) is bandlimited to W Hz and the switching 
occurs at times U =- nT + *. where T £ 2W and <„ is small com- 
pared to 7\ 

(it) The filtering aspect of the filter is neglected in the sense that if 
no switching is performed, then v(t) = x(t). In particular this 
means that if current impulses of area x(nT) are applied at times 
nT to the listeners capacitor, x(t) will occur at the output 

(iii) Suppose that one volt is placed on the listening capacitor when 
no energy is stored in the filter. The voltage across the capacitor 
(t £ 0) under these conditions will be called z(t) r and we assume 
f(0+) - i- 

Our goal will be to determine, for given spectrum of the input, and 
for given second-order statistic* of the <„ sequence, the output S/N. 
The proper design of zlt) will be of major concern in the analysis. 

We shall call the problem just describe*! two-wire switching. 

Ill- FUNDAMENTAL EQUATIONS 

In this section an exact equation which describes the two-wire 
talker-listener switching situation will be derived. As in Fig. 1. let v(t) 
be the actual voltage on the talker's capacitor and w{t) the actual 
voltage on the listener's capacitor. In the absence of switching, v(t) = 
s(t). In addition, at times nT + e. voltages \xr{nT + t n -) - v(nT + 
<b _)] are placed on the capacitor, where u(t-) means the limiting 



•Of course in reality they are not interchanged. What happens id that at the 
given instant, of time at which the "switch" is to occur, the instantaneous voltages 
are measured, certain current sources (not shown) are activated and a fixed value 
of current flows for a duration just sufficient to effect the interehnnae of chow*, 
and hence voltage*, of the equal capacitors. All this occurs in a negligible period 
of time compared with the response time of the filten*. 
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value of w(f) as t* approaches t from below. Since switching occurs at 
times nT + «„ t continuity cannot h<- assumed at these instants and left 
and right limits must be distinguished. Uaing the definition of z{t) we 
have by the superposition principle 

v(f) - x(t) + Z [w(nT + «„-) - v(nT + u -)Ht - nT - u) (1) 

where l(t) is the integer which satisfies 

KOT + ti £ t < (i(t) + \)T + *,.». 
Likewise for w(t) we may write 

If II 



W 



(J) = £ HnT + «.-) - w(nT + ..-)&</ - nT - •.)■ (2) 



If we let t — ► (kT + €*— ) in equations (1) and (21, we have the 
pair of equations 



i'o 



v(kT + «,-) = x{kT + «.)-£ HnT + *.-) - w(nT + «„-)] 

■f[(fc-n)r + e, -e.-] (3) 

and 



Mm 



■t[(k-n)T + „-«.-]■ (4) 
It is very useful to now introduce the vectors Y, V, TF and a matrix Z 
defined by 

y, = z(kT+ tl ), 

v t = v (kT + , k -), (5) 

w t = w(kT + «.-), 
and 

Z». - zl(* - n)T + «, - .--]. (6) 

Note Z kn = if k g n since z(0 = for ( < 0. 
Using equations (5) and (6), (3) and (4> become 

V - Y -Z(V -W). (7) 

W -Z{V - W). (8) 
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Solving the above pair for (V - W) gives* 

(7- W) - [I + 2Zy l Y. (9) 

Equation (9) determines exactly the behavior of the switch, for the 
output is determined by applying the sequence {\\ - W k ) at times 
kT + Q to the filter (here assumed ideal! . We emphasize that the 
jitter enters not only through the time at which the ^-functions are 
applied, but also in the quantities Z and Y. 

Consider now a z{t) such as that shown in Fig, 2 which passes 
through zero at all times kT, k > 0. If further there is no jitter, then 
Z = 0. and equations (5) and (91 show that impulses of area x{t\T) 
are applied at times w7\ thus giving x(t) at the output. In short, in 
the absence of jitter, any z(t\ which passes through zero at all positive 
integer multiples of the sampling interval will be optimum. One might 
now conclude that a zU) which is zero in a sufficiently large neighbor- 
hood of each such crossing will not see the jitter and would therefore 
be optimum in the presence of jitter. However the following argument 
(by J, F. O'Neill) suggests that such is not the case. Consider sampling 
a dc signal at random times. We must get enough power through the 
switch to reproduce the signal. Surely if we sample late, we are lagging 
in power and we would like to increase the area of the impulse; like- 
wise if wc cample early, we seem to be supplying extra power and so 
should decrease the impulse area. There should be a design of z{t) 
which would conspire with the jitter to reduce jitter noise even below* 
that noise obtained for the class of z(() which do not see the jitter, 

IV. WHAT Ifc THE OUTPUT NOISE? 

Define the vector 

y m V - W (10) 

and the functions 

sin£(/-nT) 

*-(0 = — - (U) 

f « " nT) 

for all integer n. The properties of $AD that we use are 

WO- <M-0. (»2) 



* The inverse of (/ + 2Z) always exiata since it i* triangular and all its diagonal 
element* are unity. 
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Fig. 2— An illustrative curve for *((>. 



/: 



WfliUO dt = Ti HW , , 



and 



Mi - « = E #..-('.)*-<(»)■ 



(13) 



(14) 



We assume that the reconstructing filter impulse response is given by 
$o (t) and thus the output error signal is 



and has average power N where 



(16) 



(I«) 



In the right side of equation (16), we have introduced e k (t) which 
is the error signal truncated to K pulses; that is, the upper limit of the 
sums in equation (151 is K instead of infinity. If we use equation (14) 
to expand $ m (t — e B ) which occurs in equation (15), find e*{l) and 
do the time integral, we obtain 



/_>« 



T E t-7i E fc40fc+_«fo) 



* . * ""■ 



(17) 



*tt-t 



•This corresponds to the noise in a band up to 1/2T Hz. If T is Ims than the 
Nyqiiist rate, then some out-of-hnnd noi*e is included here. In practice T will be 
about half the Xyqui&t interval in order to make filter design problems easier, 
but then the inclusion of the out-of-band noise seems fair at this stage. Figure 3 
shows a picture of the relevant bandwidths. 
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Equation (14) may be used a&ain to rccombinc the first terms of the 
right side of equation ( 17* and so 



KTJ- 



■ JC i * 






Recall again that 



Z.. = *!(* - n)T + «, - gj, 
7l = (1 + 2Z)«Vi , 
y, - x(*T + «,). 
It is useful to define a new matrix $ by 

2Z 



(19) 



that 



1 + 2Z 



«r. 



(20) 



(21) 

Substitution of equation (211 in equation (181 splits the expressions 
into two types of terms. The fire! type, called collectively A K . do not 
involve the filter z(i) while the remaining ones, called B K . do. We 
thus have 

N = lira filial + lim E[B K \ - A + B, (22) 

where in equation (22), the expectation is taken over both the signal 
and jitter statistics. The quantities A t and B t arc given by 

Ak - js £ Mft*.-.<«. - O + £ E «: - | £ I.B.UO. (23) 
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Fig. 3 — Relevant bandwidlhs in terms of the sampling interval T. 
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2 K 

+ 17 £ (eyU9y).4>*->(<> - «■)• (24) 

One may note that if »(*) is such that when the filter is designed 
so as not to feel the jitter (as described in Section III), then B = 
and the noise would be given by the filter independent torm A alone. 
The function of the filter design is to make 6 as negative as possible. 

To proceed further we make essential use of the fact that the * t are 
small. In addition* since z{kT) = is an optimum solution in the 
absence of jitter we shall keep the requirement 

zikT) = 0, *- 1,2, ... (25) 

and see what further optimization can now be made. This will amount 
to designing the slopes of the function z{t) when it goes through zero. 
The evaluation of A and B for small t is carried out in Appendix B. 
Introducing the correlation function of the signal 

R(t) = E{x(t)t(t + r)J (20) 

and the correlation function of the jitter 

J(k - n) - £[<*«.] (27) 

we find that when the jitter is independent from sample to sample, 

A = lim E[A K ] = aj[-ft(0) + | ^J?<0»]- (28) 

In equation (28), g\ ia the variance of the jitter and is given, for inde- 
pendent jitter, by 

,J = 7 - i 2 = 7(0) - J(n * 0). (29) 

To write the corresponding expression for B k we introduce, as is done 
in Appendix B, the derivatives i, of the function g(t) at the zeros, 
that is 






• - 1.2,3,.-. (30) 
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and also the constants* 

r. = ^.(0) = *«,(-*) = ^^- , a * o. (3D 

Then, for jitter independent from sample to sample, 
B - lim E\B K \ 

= 4ff(0)<rJ 



(32) 






i> c 



Expressions for .4* and #* which include the effects of correlation in 
the jitter are also derived in Appendix B and will he discussed in 
Section V- At tlie moment we merely state that positive correlation 
l>etween jitter values will tend to reduce the noise power given by the 
sum of equations (28) and 132). Equation {32> Rives, incidentally, 
more of a physical interpretation as to the breakup of the noise into 
-4 and B terms. The filter mentioned at the end of Section III which 
does not see the jitter certainly has i(t) = at each crossing, and thus 
from equation (321, B = for that filter. Thus the A term is the 
noise power for the "blind" filter; it will be improved upon whenever 
B is negative. 

If we now define the functional F\z\ of t(t) by 

TO = ir,. + 2 p. - ±%i. + 2 &*%*. (33) 

i>J 

then optimum choice of any it is found by simple differentiation of 
equation (331 : 

*i. -°- r ' + 2*- aro + 3 £f* k(0) ' *" 1 ' a - • 

(34) 

Thus for given signal statistics, equation (34) must be solved for the 
optimum set of (*?*}* In reality, equation (34) does not have to be 
taken seriously for all positive integer k, since the response of a 
realistic filter will die off rapidly with time. Thus in equation (331, 
most of the i* can be set to zero and only a few retained. If the first Z 
are thus retained we shall refer to this as designing the first Z zero 



•Recall the definition of iA*<0 given in equation (11). Also in equation (31), 
dot* denote differentiation with rc&pect to the time variable* 
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crossings* Of course equation (34) will hold only for the k such that 
iic¥=Q<x priori. 

As our first example we consider the case where the signal spectrum 
is flat and extends to Qw = w/T rad/sec. Wc refer to this as the 
full spectrum case. The interval T in this case is also the Nyquist 
interval, and no ovcrsampling is done as one would expect to do in 
practice (the latter case will be treated shortly). We have 

Full Spectrum Case: 

R(s) = 0, 9 * 0, 

fi(0) = 0, (35) 

fi(s) = B(0) ^^f , a * 0, 

Thus equation (34) yields as the optimum solution for designing all 
zeros in the full spectrum case 

4 " m ■ 2 * (36) 

Proceeding to evaluate A and B, we have 

A -fl(0)-||iff!, 

37 (37) 



b- -«(o)|p^. 



Since R{0) is the signal power, we have 

An important fact can be gleaned from equation (37) : one should not 
generally expect the optimum filter to be very much better than the 
"blind" filter/ In the reasonable example which we have just worked, 
only a gain of 3 dB is achieved. 



♦Of course the blind filter is itself & highly designed filter. 
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Table I 


—Full Spectrum 


Output S/N 


Jitter Standard Deviation 
0* 


15 dB 
21 dB 
29 dB 
35 dB 
55 dB 


0.1 T 
0.06 T 
0.02 T 
01 T 

o.ooi r 



Table I calculates equation (38t for several values of jitter. For an 
output S/N of 30 dB, two percent jitter is required. 

V. EFFECTS OF SIGNAL CORRELATION 

Equation (34), the basic design equation for the filter response 2(0, 
can be conveniently rewritten in matrix notation 

(/ + M)i = -4 (39) 

where the vectors £ and A have components 



Ak " 2 U(0) r \| 



(40a) 
fc - 1.2, ... . 



and the matrix M is defined by 



M., = S|j_2. (40b) 

The matrix / is of course the identity. The dimension of all the above 
quantities is Z, where Z is the number of zero crossings that one wishes 
to design for. 

To solve equation (391 one must in general invert the matrix 
(/ + M). This was possible in the full spectrum case because M was a 
diagonal matrix; in general it is hard to do exactly, unless the dimen- 
sion Z is small. A case of special interest for applications is a qualita- 
tive understanding of the situation for a flat signal bandwidth up to 
(Inn = ir/27\ This correspond? to sampling at twice the Nyquist rate 
and better approximates the situation to be encountered in practice. 
We call it the half spectrum case. We have 



9U4 
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Half Spectrum Case: 



ft(*T) - ft(0)^sin|s. a*0; 



/far) = rt(0> 






2 sin 2 * 



The /! term yields, from equation (28) 



(41) 



m = i2T i ° 1 ' = T' (4 ' m - 



(42) 



In the Ii term wc must solve equation (39), We do so by inverting the 
malrix by hand foi 1 the cases of designing either the first zero, the first 
two zeros, or the first three zeros (Z = 1. 2, 3 respectively). We obtain 



-M = ±1-0.409]. 



up l 









-0.454 

+0.I45J 

-0.442 

+0.179 

-0.118 



These numbers give for B 



W), = :,f i( ~ lM) 



a 



^(-1.49). 



«(0) ,-, T 
H = S(-1.5Q>. 



(43) 



(44) 



«(0) ,-, 7" 

From equations (44) and (42) we note that Z = 1 design improves the 
"blind" filter by only 17 dB and Z = 2, 3 by 1.8 dB. Output S/N are 
given in Tabic II. For a given fraction of jitter essentially the same out- 
put S/N is obtained as for the full spectrum case shown in Fig. 1. Let 
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Table II- 


-H 


\lf Spectrum* 


Output 8 


,'N 




"jr 


Z 




1 


2 


3 




15.5 dB 








0.1 


ai.sdS 


S 




8 


0.05 


29.5 dB 


A 




A 


0.02 


35.5 dB 


M 




M 


0.01 


55.5 dB 


E 




E 


0.001 



* Th« noi#e is here measured in a bandwidth twice the signal bandwidth. 

us emphasize here that the jitter is listed as a fraction of the nominal 
sampling time and not in absolute units. Thus for a fixed signal spec* 
trum and fixed jitter in seconds, one percent jitter for the full spectrum 
case would correspond to two percent for half spectrum situation. 

The effect of positive correlation can be seen by comparing equation 
(43) with equation (36). Positive correlation tends to flatten the slopes 
at the zero crossings somewhat. Pursuing positive signal correlation to 
the utmost, we consider one more case, the case of a dc signal* In this 
case $(r) = 0, and we want to solve 

(I + M)i = -*r (45) 

where 



/ + M 



(46) 



111 

1 2 I 1 -■ 
112 1- 

Let us first consider the (Z x Z) version of equation (45). One may 
verify that for dc 



(/ + J#)"' = / - 



1 



U 



(47) 



Z+ I 

where U is the Z x Z matrix that has all its elements equal to unity. 
If we let V equal the Z dimensional vector which has all its com- 
ponents equal to unity, then we get 



4 __r . _! ff (-D-'1 r 

2 T 2(Z+ I)T Ln? n . 



(48) 



(49) 
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Finally for dc signal we have 

R(Q) " T 5 ' 3 
and 

m - *M- L 2 

which has the limiting value (Z -» cc) 

Again we see for an exact solution that only 3 (IB is gained by pursu- 
ing optimum design past the "blind" filter design. The output (S/N> 
is 3 dB better than equation (381 ; one would expect dc to be influenced 
less by jitter. 

We note the remarkable fact that as Z -* co, equation (48) yields 
exactly the same solution k k = -r fr /2, as one gets for the full spectrum 
case. This is not true for a fixed finite Z. 

Judging from our two exact and one approximate solution, designing 

* = -£ (51) 

is an optimum design independent of signal spectrum. Also the "blind" 
filter is a very good design independent of signal spectrum, being about 
only S dB worse than optimum. This all assumes that the jitter is in- 
dependent from sample to sample. 

VI, EFFECTS OF JITTER CORRELATION 

Until now, any correlation between jitter samples has been ignored. 
In fact, some positive correlation is to be expected in the jitter statis- 
tics. We do not feel it is large for the way we have described the jitter 
in the previous sections* (call it *-jitter)> so that our previous design 
is not affected. Nevertheless, we should note that any positive correla- 
tion which exist* between the jitter values will improve the output 
S/N above that obtained by assuming that the jitter from sample to 
sample is independent. The physical argument as to why this should 
be true is quite simple. Assume *i is a nonzero constant, the same con- 



*8ee Section VII. 
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slant for all i\ Then we are sampling at time t n = nT + <. From equa- 
tion (19), Z = and we are supplying current impulses of area 
x{nr + c) at times nr + c to the output filter and the signal is recon- 
structed perfectly. Positive correlation helps. 

An important check on our work thus far can he made by assuming 
t - const., and checking that the -1 and B terms, which for jitter are 
given in Appendix B, vanish* Since for this case 

J(0) - J(s) t for all «, (52) 

the B term given in equation (112) obviously vanishes. The explicit 
demonstration of the vanishing of the A term in equation (107) is not 
at all so obvious. Term (107), the A term, for constant jitter, is seen to 
be proportional to (note ^ (0) ■ 0] 

-#(0) - EftttM + 2 Z*(*)iM0) 

= -*«*> - E m)M*T) - 2 E*0*o(«T). (53) 

* • 

Define a function of time g(t) by 

fft- Jim*!*). 

Then the right side of equation (53) is 

and by Poisson's sum formula 

E 9(*T) ~ \, _£ g(£ «) (54) 

where G'M is the Fourier Transform of git). Now $ (t) is band- 
limited to H < ir/? 1 and so is R(t). The spectrum of R(t)+At) ex- 
tends only to 2w/T and in fact vanishes at *> = 2w/T since it is a con- 
volution. Thus only the m = term of equation (54) could possibly 
contribute, but this contribution also vanishes because the second time 
derivative in the definition of git) introduces a double zero in the 
spectrum of 0(#) at •# = 0. 

Further use will he made of the expression for output noise with cor- 
related jitter when the «-modcl for jitter is compared with another 
model discussed later in Section IX. 
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VII- TIMING STABILIZATION 

At this point we return to consideration of many talker-listener 
pairs. Isct L be the number of possible such pairs. This number depends 
on how rapidly each talker can be processed (switched)- For example, 
if each talker is sampled at a nominal rate of \/T and t t is the process- 
ing time for the tth speaker, a constraint which reads something like 



£*. «r 



(55) 



*-i 



must be imposed. Further, the timing errors << are hopefully small, A 
scheme for stabilizing the sampling rate has, in fact, been proposed 
but T. G. Lewis.*' 2 An interpretation of this scheme due to Saltzberg 
and Pasternak is shown in Fig. 4, For immediate convenience let us not 
worry about normalization and let the dashed line have unity slope. 
The large dots in Fig. 4 represent a talker starting to be sampled, and 
let the length of time required for the ith sample be h . The time U 
represents the length of time that the current source pumps current to 
effect the interchange of the talker-listener capacitors and, in accord- 
ance with our earlier assumptions about measurement times being 
short compared to filter bandwidths, we have * 4 <K 7\ After many 
counts (measurements) , say L r we return to resample a given talker. 
We particularly note that a talker is never sampled early. We see that 
in Fig. 4 after the time t s , our path hit* the 45° line early, and we 
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Fig. 4 — Timing stabilisation. 
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K p»d" time until the fourth count occurs. Also, as shown, the variables 
q n represent the horizontal distance from the 45° line to the dot repre- 
senting the count. They are positive random variables and are given by 
the equation* 

q 9 - max J (56) 

Ifr^ + t-i - 1. 

The random length of times t*-\ are assumed to be identically dis- 
tributed and independent of q n , x . We shall proceed to derive an in- 
tegral equation for the probability distribution function for q H S We 
note the possibility of a S-function at q m = and thus write for the 
density function g n (x) 

q m (x) - *. «(x) + (1 - cOp B (x). (57) 

The number a, is the probability of having q n - 0, and p,(z) describes 
the continuous portion of the density. If we let u(t) denote the density 
of t* and * denote convolution, we have, using equation (56), 

a n = *,-, f u(t) dt + (1 - a,-,) | [p.., * w] dt\ (58a) 



*» = «-. ^ jtJL -(!-«.-.) 



f~ p^{t*)u(x + 1 - t')dt' 

f u(t') dt' f~ dt j~ p.-i(*>(/ - /') dt' 

(58b) 

A steady-state a and p{x) would obey equation (58) with all indices 
removed. We shall write the steady-state equation. Let K denote the 
known constant, 

K - fu{t')dt' f (59) 

let, for x > 0, v(t) denote the known density function 

,(0 = U -^± . (60) 



•Except for a time scale normaliiation. the sequence {q mL } n = L 9 ■■■ conr- 

Bonds to the sequence M of previous sections when there are L talkera on 
e switch. 

* B R. Salteberg and G. P. Pasternak 3 had also derived an integral equation 
for ifc, and no doubt our equation it aubsunnally equivalent to their equation. 
An entirely different approach (i*>, an equation for a different set of variables to 
describe the same problem) haa also been discussed by J. Salx and R. D. Gitlin.* 



■ - *ty m 

One may now imagine determining the unknown constant P in the 
following manner. Equation (63) determines CM and therefore tjie 
density p(x) in terms of p. Form the convolution of p with u and set 
the finite integral of this convolution equal to p in accordance with 
equation (61). This then is an equation which can be used for the 
numerical determination of p. For the case 

«(0 - exp (-fit) <65) 

the procedure may be carried through exactly. We find that the dens- 
ity p(x) of the q variable is also exponential and is 

p(x) - <*0exp(-a0x) t (66) 

where the probability a of having q exactly zero is related to by 

1 - a = exp {-0a). (67) 

Using the inequality 

1 - x < e'*, x * 0, 

we see that equation (67) has a nonvanishing solution for a if and 
only i! > 1. Realizing that the average of /< is 1//?, this says the 
system will be stable if the average duration of U is less than one. 
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and introduce the unknown constant 

P= f \p*u(n)dt'. (61) 

Jo 

Then the steady-state solution (lira n -* oo) satisfies, if it exists, 

a= (1 -K)a+ (1 - a)p, (62a) 

p(x) - av(x) + (1 - a) fp(On(x - t') dt\ (62b) 

Jo 

If we let CM denote the characteristic function of p{x) then equation 
(62) yields 

«« - . - (fJ"W) (98) 

where 
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From the above consideration?* we calculate that 

VIII. EFFICIENCY OF SWITCH 

We propose here to give an idea of how much is gained by asyn- 
chronous (but stabilised) versus synchronous sampling. We have 
shown in the previous section that (introducing unnormalized quanti- 
ties now) for a speech processing duration distributed exponentially, 
i.e M 

«(t) = fie'", (69) 

we have 

. 1-J 



*. ■ — «*- (70) 

a p 

where a and /? are related by 

I - a = exp [-Sf£] . (71) 

which requires 

8 > \< (72) 

Again, h is the number of talkers. To effect the comparison with the 
synchronous version of the above scheme a new parameter has to be 
introduced which represents the peak-to-average voltage (not power) 
ratio of the signal. Recall 1/0 is the average processing time, and 
therefore represents the average voltage. In the synchronous case a 
maximum time, t mnx , is allowed for each talker. Clearly In, is a 
representative of the peak voltage. The ratio -4 is then 

A m 0U, . (73) 

If /,, is the number of synchronous talkers, we also have 

L.U. = F. (74) 

Now u?e equation (70) to solve for a, substitute the expression for a 
into equation (71), eliminate fi via equations (73) and (74) and obtain 



1002 THE BELL SYSTEM TECHNICAL JOURNAL, MARCH 1971 



1 . 4 « exp 



-fe) 



(75) 



where 

An approximate solution to equation (75), which is accurate for large 
L B , is that 

r-**A, (77) 

i.e., in the large L& limit, the ratio of possible asynchronous talkers to 
synchronous ones is approximately the peak-to-average voltage ratio 
of the speech signal. One would expect this number to be at least 
four. Actually for L& of interest (77) is not an accurate enough ap- 
proximation. The solution of equations (75) and (76) is shown if 
Figs. 5 and 6 for A =4 and 8 respectively. Taking A = 4 (which i- 
conservative), wc see from Fig. 5 that if the technology would permit 
50 talkers to be put on the switch synchronously, the switch could 
accommodate 150 simultaneous speakers, sampled on a stabilized asyn- 
chronous manner, with an output S/N of 35 dB.* To repeat, factor of 
3 increase in efficiency seems like a conservative estimate. 

IX. FOUR-WIRE CONSIDERATIONS 

The last topic we consider is what we call a four-wire treatment of 
the problem. Here we treat the details of a model proposed by F. K. 
Becker* which should deal with active asynchronous energy transfer 
when four-wire facilities are available. The model is this. A speech 
waveform is sampled and held for a variable time A< before the next 
sample is taken (see Fig. 7). The holding times A< arc independent, 
identically distributed, and have average value 

A. = T. (78) 

The initial speech waveform is approximately reconstructed by passing 
the jittered box car through a filter having an inband characteristic 

A(*)--Vx P (f). < 79 > 

siny 



•The 9/N ia read from Table II for *,/T - 001. 
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Fig* &— Efficiency versus number of synchronous talkers for peak-to-uvrrage 
voltage ratio 4=4. 



The Af in the above model will not necessarily be assumed to be 
sharply distributed about T\ later a stabilized version using an t-jitter 

model will be used- 
Let x(t) denote the speech wave and Xjit) denote the jittered box 
car version. Further denote the spectra of the two processes by S(«) 
and Sj{*>) respectively. Our major interest shall focus on determining 
the spectrum E(<*) of the error 



error = h * x,(t) - x(t) t 



(80) 
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Fig. &— Efficiency versus number of synchronous talkers for peak-to-averagc 
voltage ratio A = 8, 
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t t, t* t a 

Fig. 7— "Jittered-boxcar" sampled speech wave. 

assuming h(t) is the reconstructing filter. The output noise power can 
then be obtained by integrating the spectrum. It may be shown that for 
the d-modcl of jitter described above the Mart error spectrum is 



1 - ^ **<»>[!- C<*>] 



E(») - £(«) 

+ HF»{ P -e W £™ T ±^g}.« 

In equation (81 ) , the function C U) is the characteristic function of the 
variable A, that is, 

C(o>) = [' cxp (wA)p(A) dA, (82) 

p(A) being the probability density of 4. Also 

A(«) - f MP (-toO*<0 <**• 

The symbol P in front of the integral in equation (81) denotes that 
the principal value is to be taken for the simple pole 1/| 1 — €[+>' — 
to)]. If the term in the braces in equation (81) is rewritten as 



-L 



S(u') Re 



2t 



=-, (83) 



(1 - C*fo))(l ~ W)) 
. 1 - C(w' - u) . 

the reader should be assured that now no singularity will arise in the 
integrand, and the principal value distinction need not be made. 

For numerical purposes, we plot the error spectrum given in equa- 
tions (81) and (83) for the case when A is gaussian distributed about 
mean value T, and for standard deviations (*&/T) = 0.310. 0.1.* 



•The random variable A is always positive, while the RAuesian assumption al- 
low* it to become negative with some probability. We have verified that this 
effect ia oot significant here. 
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These larger variances are chosen here because we are mainly in- 
terested in the unstabilized version. Figures 8 and 9 show, respec- 
tively, these error spectra, assuming the input spectrum is flat up to a 
maximum frequency n = t rad/s, and assuming a sampling interval 
r = 05 8, This corresponds to sampling at twice the Nyquist rate. 
Also in calculating the out-of-band noise we have taken the reconstruct* 
ing filter (79) to extend to twice the baseband spectrum. This em- 
phasizes the higher frequencies more than desirable, perhaps* Thus 
in Fig. 9, the S/N for noise measured in twice the baseband interval is 
24 dB while inband we have a S/N of 29 dB. Thus for independent A 4 
a 10 percent jitter about the mean value seems tolerable for sampling 

at about twice the Nyquist rate/ 

In the event that the unstabilized four-wire version is unsatisfactory, 
we recalculate the error spectrum when timing is stabilized according 
to c-model discussed earlier. In this case the Fourier transform e L M 
of the truncated error signal is found to be given by 

X I exp ("^ - W,) " e*p ("y - <«.) J 



*£ _L_ fofl'Vl . (^ me terms\ 



2 —JT exp VT/ I" With t = 0/ 



sin -g 



j 



(84) 



Again, the reconstructing filter (79) has been assumed. Further details 
of the calculation will not be recorded, but we simply state that the 
error spectrum EM is, retaining only second-order terms as in pre- 
vious work involving c-jitter, 







l)+R(s - 1) 



•1 



■]■ 



- sin* ^ R(s)J(s) 



•Let us remark that the A model for jitter described here is significantly dif- 
ferent from the t-model used previously. The A model implies positive correlation 
between adjacent * variables and for a Riven variance produces leas noise. We 
shall return to theae quest-ion shortly. " 't 
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|sin y J(*)[(-R(«) + £(« + 1)) exp 

¥)]}■ 



- (-R(s) + R(s - 1)) exp 



(85) 



One may verify that the above expression is indeed real, that is, it 
vanishes when ■* = t and when the correlation function corresponds 
to a dc input signal. The noise power N in any bandwidth a is gotten 
from equation (85) by 



N = ;J- f E(u) dw. 
aw J -q 

Choosing fi s */T we obtain from equation (86) 

N = -*«V(0) + ^ ± a. y W [«« - *<' + D + «f - '> ] 

I 



(86) 



T.fr 



-|;£j. '«[-** + afe ^ 



i ) + Ms - i) 



(87) 



0.04 r- 




Fig* 8— Error spectrum for T = 1/2, * A /T = 0316, Undistorted spectrum is one 
for — » ^ * ^ *. and sero otherwise, 8/Nw = 185 dB, S/N* = 14 dB. 
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Fig. 9— Error spectrum for T = 1/2. *&/T = 0.1. Undistorted spectrum ift one 
for -ff<tfSf, and »ro otherwise. S/Nr = 29 dB, S/N fc = 24 dB. 



when 



1 f* cos ax i , 



' 2 



r rfr, 



8111 2 
*, = - f* xzmaxdx - ( ~^' 

IT /„ « 



H 0. 



(88b) 



(88c) 



Numerically 

a„ = 20 o w 2.78, 
«, = -0.469, 

0„ - 1.39, 
0, = 0.386. 
For the case of inde|>endent (-jitter, equation (87) reduces to 



(89) 



noise power; 
independent 



i von T * <0 > x ! « d ?i!lUlfl 4lll 



(90) 
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The results for the stabilized full spectrum And half spectrum cases 
arc obtained by using equations (35) and (41) in equation (90). 



Full Spectrum: 

In) = r |j~ + 

Half Spectrum: 



^(4.67). (91) 



(0-*K-4-D-*JW 



(0.946). (92) 



By eomparison with equation (38), the four-wire full spectrum answer 
(91) is about 1.5 dB worse than the optimum two-wire result, but 
the four-wire result has about a 4.5-dB advantage for the half-spec- 
trum case [compare equation (92) with equations (42) and (44)]. 
The latter case corresponds more to the case of practical interest. We 
note further that in this comparison the four-wire might be penalized 
by two dB or so due to the (xf sin x) characteristic being used out of 
band also. Secondly, the more concentrated toward dc the signal spec- 
trum is, the better the four-wire version will become since it has 
vanishing distortion for dc, while the optimum two-wire version does 
not 

Before leaving this topic, a comment on the relation between 4-jitter 
and e-jitter should be made. We begin by looking at two successive 
A-variables A„ and A n -i in terms of ^-variables. 

A.-r + *-^. (93) 

Clearly 

B(t) - T, (94a) 

*\ - 2?(1 - p), (94b) 

E[{^ x - 30(4. - T)] m 2?p ~t-7. (94c) 

In writing equation (94 a through c), we define 

£[«.«..,] = ?P (95) 

and assume no ^correlation after one displacement, i.e., 

B( Wt ) = i\ 
To make two successive a variables uncorrected, we set equation 
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(94c) equal to zero and obtain 
and 



(97) 



ff A = <r« 



Thus to compare independent A-jitter with <-jitter f a good way to 
do it is to choose the same variance according to equation (97) and 
introduce a correlation between adjacent *'s by equation (96). This 
simple way of comparing the two schemes is not exact, but should be 
good if signal correlations do not persist for extended periods. In any 
event, it illustrates why for a] = v* the independent A-model will yield 
appreciably better results than the independent e model. The first 
implies positive correlation in the second. 

X. CONCLUSIONS 

We began our study with consideration of the two-wire switching 
problem (Fig, 1). After obtaining a concise mathematical description 
of the operation of this switch, (equation (9)] f attention focused on 
the optimum capacitor discharge z(t) when one volt is placed on the 
capacitor (Fig. 2). We have seen that a good design for the slopes i k 
of the function zit\ as it passes through its zeros (see Fig. 2) are the 
values i k = (-l)V(2fcT). Exact solutions for full spectrum and for 
dc indicate that this result is not sensitive to the signal spectrum. This 
optimum design can be expected to yield only a 3 dB improvement 
over the "blind" filter which if defined to be flat at the zero crossings 
and consequently docs not "see" the jitter. Typically two percent jitter 
yields an output S/N of 30 dB. Any positive correlation in the jitter 
will tend to improve this figure- 
When the processing time of an individual speaker is exponentially 
distribute*!, an exact distribution of the timing jitter with a stabilized 
clock is obtained. This result is used to study the efficiency of the 
asynchronous switch. This increase in capacity over a switch using the 
same technology but employing a synchronous strategy is, roughly, 
the peak-to-avernge voltage ratio of the signal. More accurate descrip- 
tions are presented in Figs. 5 and 6. Conservatively, the asynchronous 
switch should handle three times the traffic- 
Finally, a four-wire version modeled as a jittered boxcar recon- 
structed with a filter has been considered. Under the independent 
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4-modeI for jitter, Fig, 7, 10 percent jitter can yield 30 dB S/N, largely 
because of the correlations this implies for €-jitter. Exact error spectra 
have been plotted for this c -.-■■ also (Figs. 8 and 9). With timing 
stabilisation, and consequently the c-jitter model, the reconstructed 
boxcar should yield several dB improvement over the two-wire results 
for cases of practical interest (spectrum concentrated at dc). 
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APPENDIX A 

Some Lemmas on Limits of Sums 

We list here certain lemmas which will be needed in Appendix B. 

Isemma I: If either ^7-o I /(*) I or S*o »/(•) converges,* then 

lim| fe /(n-k) = £/(•). (98) 



n£* 



Lemma Z: If X?-o /(*) converge*, then 



Likewise 



Umi E fe«»>- fe/«. (99) 



I'm ^ E E /(». = E /(». 0- (100) 



Lemma 8: If £ /(j\ *) concern, tften 






"mi E E E /0*.«) = E /(j,*>- (101) 



APPENDIX B 



Derivation o/ Owfpu* iVoise Pouter 
Our first step is to derive the filter independent noise, i.e. 
* Neither condition implies the other; each implies the convergence of £/(«), 
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A = lim A K 

where A K is given by equation (23). Our expressions are valid to 
second order in the {<«}, Wc note the following relations 

y* = x(nT + O fc x(nT) + * n i(nT) + i&faT); (102) 

» ^(Q)i, + IJJDtfA if m ^ 0; (103) 

^(0) = -iMO); 
making use of this in equation (23), we have, to second order in {e ( } 

K -• (104) 

" a. 4-0 

Averaging over the signal statistics by using equation (26) and the 
further relations 

B[t(t + r)z(t)] - Mr) - -E[±(t)x(t + r)L ^ 

g[J(( + r)i(f)] = -#(r), 
and also averaging over the jitter variables using equation (27) gives: 



ea k - ^ z [-&mm - Mowmw)) 

-J I \R(n - k)J(n - m n -M - 2r._»J(n - k)R(n - ft)]. 

Observing that the second sum is a symmetric function of {n — k) 
and using Lemma 1 of Appendix A, the required limiting operation 
yields 

A - -ft(0)J(0) - £ «WJ(«)W0) + 2 E r.J(«)A<«). (107) 
We list for further possible 



mi 

use 
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MO) = 2 £=JjP , t fi 0; 

8 1 (108) 



&(0>« -|£; 



and a|§o 



Z #.«D " 0. (109) 



The latter relation may be derived by direct evaluation of the sum, of 
course, or by using the Poisson sum formula and the bandlimited 
nature of ^(J). 

We now proceed to derive an expression for the filter dependent 
terms B, starting from equation (24). Gathering the second-order 
terms proceeds much as it did for the A k terms, with the following 
additional complication. Recall from equation (20) that is given by 

= 2£- 42* + 82*- ■•■ 



1 + 2Z 
From equation (6), 

= 0, if k < n. 

Thus the matrix is at least linear in the {«,). Proceeding now to 
collect terms and averaging, we have to second order 

K[Bk] -IE t L-A-,R(l - m) 

[J(0) + J(m - - J(n - m) - J{n - Q] 

- I £ Ei.-fi(n - m)[./{0) - y(n - m)] (ll0) 

+ V S Zfr„_A_ 1 .ft(n - m)[J(n - A) - /(« - «)]. 

In the second term of the above, we change summation variables by 
replacing the sum over m with a sum over $, where 

a = n — m, « — l t 2, --- i ft. 
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Wo make the same change in the first term, and also introduce 

t = n - t, t = 1, 2, ■•■ , n. 
In the last term let 

& = k — m f $ = 1,2, * - - t k t 

j = n-k, j=-<*- 1), •■■ ,K -k, 

Then 

E[B K ] = £ E E i.i,R(« - OtAO) + /(* - - /« - A01 

-4e t,*M)[m-m] a") 

+ 1 E E E r^ffo + »)[J(j) - J(i + *))■ 

Making use of Lemmas 2 and 3 of Appendix A yields 



B = lim E[B K ) = -4 £i,ff(s)[J(0) - J{s)J 

+ 4 E «A*<« - 0[^(0) + J{t - t) - J(a) - J(Q\ (112) 

+ 4 E E r,i.R(j + s)[j(j) - J(j + *)]. 



#-i »--* 



The simpler equations (28) and (32) which hold for independent 
jitter may be Rotten from the above by using the result of Section 
VI. Namely for J(s) = const., all &, the -4 and B terms vanish. When 
we have independent jitter J(s) = const,, s # 0. So we add and sub- 
tract a term with J(0) replaced by J(s ¥- 0) and use equation (29) to 
obtain the independent jitter results. 

appendix c 

Exact Error Spectntin for Independent ^Jitter 
Consider a long time L f and denote by [L\ the best integer so that 

L&[L]T. (113) 
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The jittered box car (truncated to the long interval L) can then be 
represented as 



I 

a-0 



*/(0 = Z*(U5[UiiU, (H4) 



where B is a unit box extending from t. to t n .i . The Fourier transform 
of this truncated version is then 

s t («) - l E <Q «P (-i«(,)[exp (-w A,) - 11- (115) 

The power spectra for a process y(t) is generally calculated according 
to 

<7<«)« limf E|/G»)l\ 

where y £ M is the Fourier transform of y(t) truncated to an interval 
L f and E denotes an ensemble average with respect to all random 
parameters. If we apply the above to spectrum of the error signal 
(80) ^ we obtain immediately, using notation of the text, 

E(*>) = | fc(«) T 8A») + 8(0) + D(«) (U6) 

where 

/)(«) - - Mm TTvr^E Re i^JWM. ("?) 

Expressing x L (*) as a time integral and performing the signal averag- 
ing we get 

2 r •' lt ' 



DM = -lim ^^ Re ft»<- £ exp (io-Ulexp (i» AJ - 1] 

(118) 



Kir -i« .-o 



x/"exp(-toJ')B(«'- ««}■ 



Expressing the signal correlation function in terms of the spectrum, 
doing the integrals and remaining averages, gives 

DM - =r Re **(»)(^)s(«)fl - C(«)]. (119) 

There remains the calculation of Sj(*>). Using equation (115), we have 
(performing averages over signals) 
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Ei.-Ii. 



-#-o 



X [exp (-w A.) - IKexp (w A_) - l][exp - to(g A. - g A,)J 

(120) 



whore we note the fact 



11-1 



Ea.. 



(121) 



1-0 



Performing the A-avcraging yields 

[lit - mn? 2[1 Re C(w)1 k l 



III ((I 



+ — ^Re [1-C»1 S E 

/I & 5 <*W " *H<V " •* " C(w/)I }' (122) 

To perform the limits of the sums we first replace C*-*-'^ — ») by 
[exp(— «)C^"™" 1 («' — *n)l to insure convergence and then let * — ► 0. 
We then have 

in— lAJi 



Re 



tf 



{„-c(„) 1 /;j^) ,_v^-l (,23 > 



Making use of the general relation 



x ± 



1 U "'I^^Wi 



where P denotes principal value, finally yields 

W I 

+ ^ >Re { [1 _ c . M]P jrg: SM _L^_] 

The interested reader may indeed cheek that for 

C(w) = exp(u*T), 



(124) 
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the above expression reproduces the known result for a wave-form 
impulse modulated at precisely ^-second intervals* 
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